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Abstract 

We study the dynamics of diffusing particles in one space dimension 
with annihilation on collision and nucleation (creation of particles) with 
constant probability per unit time and length. The cases of nucleation 
of single particles and nucleation in pairs are considered. A new method 
of analysis permits exact calculation of the steady-state density and its 
time evolution in terms of the three parameters describing the microscopic 
dynamics: the nucleation rate, the initial separation of nucleated pairs, 
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and the difTusivity of a particle. For paired nucleation at sufficiently small 
initial separation the nucleation rate is proportional to the square of the 
steady-state density. For unpaired nucleation, and for paired nucleation 
at sufficiently large initial separation, the nucleation rate is proportional 
to the cube of the steady-state density. 

1 Introduction 

Reaction rates controlled by collisions between diffusing particles depend on the 
distribution of distances between particles as well as on the density of particles. 
In particular, as Noyes stated in 1961 ^^Any rigorous treatment of chemical kinet- 
ics in solution must consider concentration gradients that are established by the 
existence of the reaction itself n Here, we study the dynamics of point particles 
in one dimension, nucleated at random positions and times, then diffusing until 
colliding with and annihilating another particle. Competition between nucle- 
ation and annihilation produces a statistically steady state with a well-defined 
mean density of particles and distribution of distances between particles. We 
shall contrast two types of nucleation: unpaired, in which particles are deposited 
at random locations at random times, and paired, in which pairs of particles are 
deposited at random locations. The dynamics is as follows: 

(i) Particles are nucleated in pairs with initial separation b; 

(ii) Nucleation occurs at random times and positions with rate F; 

(iii) Once born, all particles diffuse independently with difFusivity D; and 

(iv) Particles annihilate on collision. 

A portion of a typical realization of these dynamics is shown in Fig. |^. For 
unpaired nucleation (i) and (ii) are replaced by 

(i') Particles are nucleated at random times and positions with rate Q. 

An existing method of analysis, based on a truncated hierarchy of correlation 
functions, is developed and extended in this article to the case of paired nucle- 
ation, yielding expressions for the correlation functions in the steady state, and 
for the time scales for relaxation towards the steady state. We also introduce 
a different method of analysis that yields an exact explicit expression for the 
steady-state density and for the time dependence of the density starting from 
arbitrary initial conditions. Our analytical predictions are compared with the 
results of direct numerical simulations. In the simulations, large numbers of 
diffusing particles are simultaneously evolved in continuous space, with annihi- 
lation whenever two paths cross and nucleation (paired or unpaired) at random 
times and positions. 

A striking difference between paired and unpaired nucleation is the scaling 
of the steady state density of particles, po, with the nucleation rate: po cx 
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Figure 1: One part of a numerical solution with paired nucleation, governed by 
(i)-(iv). Time increases upwards and each dot indicates the space-time position 
of a diffusing particle. T = f .25 x IQ-^, & 2, D = 0.1 



Fa (paired) versus po oc (unpaired). Here, we shall exhibit the crossover 
between these two cases in terms of the following dimensionless quantity: 



b. (1) 



For £ ^ oo, the dynamics described by (i)-(iv) is equivalent to that described 
by (i'), (iii)-(iv) with the replacement 

Q 2T. (2) 

The paper is arranged as follows. In the remainder of this section we sum- 
marize published results for reaction-diffusion systems. In Sec. |^ we analyze the 
dynamics using a hierarchy of equations for particle density functions, called 
"reduced distribution functions" by Van Kampen.l3 Derivation of the reac- 
tion kernel leads to an exact relation between the density of particles and the 
derivative of the correlation function. We also explore the linear response to a 
perturbation away from the steady state to establish the time scales for relax- 
ation. In Sec. ^, by introducing a function that satisfies a closed linear partial 
differential equation, we present exact expressions for the steady-state density 
and for the time evolution of the density with arbitrary initial conditions. In 
particular, analytical results are presented describing the rapid initial annihila- 
tion that transforms an initially random distribution into one characterized by 
an effective repulsion between particles. 
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1.1 Unpaired nucleation 

Analysis of difFusion-limited reaction dates back_to M. von Smoluchowski. His 
Mathematische Theorie der raschen Koagulatioru considered reaction between 
diffusing particles resulting in merger, with the reaction taken to occur immedi- 
ately whenever two particles are a distance R apart. He introduced a diffusion 
equation for the density of particles relative to the position of a test particle 
and noted that the density is zero at all times at radius Ra For many years it 
was assumed that the final result of a complete calculation following the proce- 
dure outlined by Smoluchowski would be an equation for the mean density of 
particles, p, of the forniB 

p^Q^hp^, (3) 

where Q is the rate (per unit length and time) of appearance of new particles 
and kg is constant. This would imply, for the case without nucleation ((5 = 0), 
that the density is proportional to t^^ for t oo. However, arguments basedjpn 
dimensional analysis and scaling show that this is not true in one dimension.Bcl 
In 1983, Torney and McConncU studied this case and published an exact solution 
for the mean density as a function of timeJlj Starting from an initial random 
distribution of particles, they found 

p{t) = p{0) cxp {8Dtp^{0)) crfc (p(0)(8L't)^) . (4) 

In particular, p — > (87r_Di)^^ for t — s- cx).^ rederivation of the result of Torney 
and McConnell was provided by Spougeji^ whose insight was that an annihila- 
tion process is equivalent to a coagulation process if coagulants made up of an 
even number of partial^ are considered as diffusing "ghosts." Derivations based 
on reflection principldlj and field theoryll3'Ej methods have also been published. 

In discrete models of diffusion-limited reaction, diffusion is approximated 
by hopping between neighboring sites on a lattice. Here, too—tlie density of 
particles without nucleation is proportional to t~2 for t ooJiinla Moreover, 
with unpaired nucleation, thcj-stpady state density is proportional to the third 
power of the nucleation ratc.tacil This can be interpreted as evidence for a 
time-dependent rate constant kg in (|^) , or as requiring (^) to be replaced by an 
equation of the form 

p^Q~k,p^ (5) 

However, no polynomial equation for the density can describe both the steady 
state .lAfith nucleation and the long-time decay of the density without nucle- 
ation.ElE3 

An exact solution has been found in one dimension for a discrete coagulation 
model with one fixed source. The latter solution is related to the probability 
that a given spin in an Ising. chain with random initial conditions does not 
change its value before time ttB For discrete and continuous coagulation models, 
exact results are available not only for the density but also for the spectrum of 
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relaxation_ijates,c3 the distribution of interparticle distances £3 and correlation 
functions-O They are obtained by considering the function E{nAx, t), defined as 
the probability that an arbitrarily chosen segment of n consecutive sites contains 
no particles, satisfying a closed kinetic equation. It has, however, not proven 
possible to extend this method to the case of annihilatioi3_on contact, because 
the function E{nAx, t) does not satisfy a closed equation^ 

1.2 Paired nucleation 

The "coefficient of recombination" of two particles initially close together was 
introduced in the study of subatomic particles.ES The relative motion of two 
diffusmg particles is equivalent to a problem of Brownian motion of one parti- 
cleBB 

A discrete model that corresponds to paired nucleation is the Ising model, 
with nucleation |at neighboring sites. Its dynamics was studied analytically by 
Glauber in 1963;E3 the nucleation rate is proportional to the square of the steady 
state density for nucleation rates sufficiently small that excluded volume effects 
can be neglected.Ej Computer simulations of a discretized reaction-diffusion 
model A + B ^ 0, published in 1987,Ej contrasted the scalings of the steady 
state density according to whether nucleation occurred at random sites or in 
pairs at neighboring sites. In the latter case, the scaling F cx was found. 

A different approach to diffusion-limited reaction was recently introduced 
in the ixintext of kink dynamics in a stochastic partial differential equation 
(PDE)E3 There, the dynamics was termed "mesoscopic" because it was an 
approximate model that ignored the internal structure of kinks and antikinks, 
treating them simply as particles that happen to be nucleated in pairs. The 
treatment was based on classifying particles according to whether they are an- 
nihilated in a collision with their nucleation partner (recombination) or with a 
different particle (nonrecombinant annihilation). The steady-state density po is 
related to the mean lifetime of a particle, t, by 

PO - 2Tt. (6) 

The mean lifetime r was estimated directly by averaging over the possible his- 
tories of a pair of particles born together. This approximate analysis yielded 
the estimate po ~ {SbT /SD)^ . 

2 Hierarchy of distribution functions 

Let fn{xi, ■ ■ ■Xn]t)dxi . . - dxn be the probability that there is one particle in 
{xi,xi + dxi), one in {x2,X2 + dx2), and oneJn (a;„,a;„ -I- dx„) at time 
t, regardless of the positions of the other particles.S The function /i(xi;t) is 
the particle density at xi at time t. On deriving the diffcix^ntiaLequation for 
its time derivative, one finds that it involves /2(a:;i, 2:2; i)-crQE3'Eil Similarly, 
the time derivative of f2{xi, X2', t) involves fsixi, X2, x^; t). One is thus led to a 
hierarchy of differential equations for the evolution of the distribution functions. 
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In this section we derive the source terms appropriate for paired nucleation 
in the hierarchy of differential equations. We also derive the reaction terms 
corresponding to diffusion with annihilation on collision, without needing to in- 
troduce a reaction radius. Three parameters remain in the theory for paired 
nucleation: the nucleation rate of pairs T, their separation at nucleation b and 
the diffusivity of a particle D. For unpaired nucleation there are two parame- 
ters: the nucleation rate Q and the diffusivity of a particle D. The annihilation 
process is immediate on collision and therefore does not require extra param- 
eters. It manifests itself instead in boundary conditions on the distribution 
functions. We shall truncate the hierarchy of distribution functions using an 
ansatz for the three-jioiiit correlation function, introduced in the literature for 
unpaired nucleationj£3H thus obtaining a closed pair of differential equations 
for the density and two-point correlation function. Their solution yields an- 
alytical approximations for the steady-state density and two-point correlation 
function. By examining perturbations away from the steady state, we derive 
the time scales for relaxation towards the steady state. 

The evolution of the reduced distribution functions has a number of contri- 
butions 

d 

/„(xi, . . .x„;i) = DW'^fn{xi,...x,i;t) 



at 



^ ^ ^(^i 1 ^j^fn (-^1 7 ■ ■ ■ •^n 7 ^) 



/ <^Xn+l k{Xi,Xn+l)f7i+l{xi,. . .Xn,Xn+l\t) 

i=l 



+ sources. (7) 

First on the right is the diffusion term, due to the motion of each particle with 
diffusion coefficient D. The second term represents the reaction between two of 
the n particles: k(x, x') is the probability per unit time that a particle at x and 
one at x' react, and the summation is over all pairs that can be selected from 
the n particles. The third term accounts for the fact that each of the n particles 
may react with another that is not part of the set of n particles, and 2L is the 
size of the system. The last term is a source contribution whose form is given 
in detail below. Equation (m) is one in an infinite hierarchy. Explicitly, the first 
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two equations in the hierarchy are: 

d 9^ ['^ 

-Q-Ji{xi;t) = D-^fi{xi;t) - J dx' k{xi,x')f2ixi,x';t) + qi{xi), 



(8) 



d 

f2{xi,X2;t) = D\/'^f2{xi,X2;t) - k{xi,X2)f2{xi,X2;t) 



dt 

L 

dx'[k(xi,x') + k{x2.,x')]f:i{xi,X2,x' -,1) 

-L 

+ q2{xi,X2) + h{xi\t)qi{x2) + fi{x2;t)qi{xi). (9) 
2.1 Source terms for paired nucleation 

The term qi{xi) in (||) and (||) is the probabihty density per unit time for the 
creation of a particle at xi] the term q2(xi,X2) is the probabihty density per 
unit time for the simultaneous creation of a particle at xi and another at X2- 
When creation of particles always occurs in pairs, these two source functions 
are related 

qi{xi) = j dx2 q2{xi,X2). (10) 

When the particle creation rates are independent of position and time, qi{x) is 
constant 

gi(x)=2r, (11) 

and q2{x,x + y) is independent of x. The constant T is the rate of creation of 
pairs per unit length. 

Here, because particles are indistinguishable, 92(2:1, X2) depends only on 
y — \xi — X2I and the functions fnixi, ... ,Xi, . . . , x„; t) are independent of the 
order of the Xi. The probability that two particles initially at xi and X2 react is 
the probability that they diffuse and collide. Since particle diffusion is isotropic 
and independent of position, k{xi^X2) also depends only on y = \xi — X2\. We 
therefore define 

q{v) = q2{xi,X2) and K(y) = k{xi,X2). (12) 

The function q{y) describes the probability density of distances between parti- 
cles nucleated simultaneously. We shall use the following forms for this function, 
corresponding to unpaired nucleation and to paired nucleation with initial sep- 
aration b: 

unpaired; 

(13) 

paired. 
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We can now rewrite (||) and (g) as follows 



^Ji{x;t) = D^h{x-t)-2j dyK{y)S2{x,x + y;t) + 2T, 

(14) 

d 

f2{x,x + y;t) = D\/'^f2{x,x + y;t)-K{y)f2{x,x + y;t) 



dt 



L 

dz {K{\z\) + K{\z - y\)) hix, x + y,x + z;t) 



+ q{y) + 2r(/i(x; t) + h{x + y; t)). (15) 

If the initial conditions are homogeneous, then the functions /i, /2, ... will 
be homogeneous at all times. In particular, fi{x; t) will be independent of x at 
every t. Let 

p{t) = h{x;t), 

Fn{x2, ■ ■ ■ ,Xn]t) = p^"^ (t) fn{x , X2 - X , . . . , Xn - X] t) . (16) 

We shall in particular be interested in the dimensionless correlation function 
defined by 

g{y,t) = F2{y;t). (17) 

The function g{y, t) is the probability density at time t of particles at a distance 
y from a reference particle, divided by the overall density of particles. It is 
constructed numerically as follows. Choose a sample of A'' reference particles, 
located at {x^, i = 1, . . . , iV} at time t. For each x^, construct 

Gi{y, t) — {number of particles between Xi and Xi + y} 

for y > 0. Then G{y,t) is the average over the N particles of the Gi{y,t) and 

g{y,t) = {p{t)r^^G{y,t). (18) 

If there is no correlation between particles at time t, then g{y,t)=l for all 
y > 0. In all the situations considered here, the total length 2L of the system 
is sufficently large compared to the correlation length so that 

lim go{y) = 1, (19) 

y— »cxD 

where go (y) denotes the steady state correlation function. 
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In terms of p(t) and g{y,t), Eqs. (14) and ( p^ now simplify to the pair of 
equations 



l-L 

„2 



p{t) = -2p\t) dy K{y)g{y,t) + 2T, (20) 







3 d'^ 

-g(y,t) = 2D—g{y,t)-K{y)g{y,t) 

- p{t) j ^dz {K{\z\) + K{\z - y\)) Fsiy, z; t) 

+ ATp-\t) + p-\t)q{y) - 2p-\t)g{y,t)p{t)- (21) 
2.2 The reaction kernel 

To complete the description of the dynamics of the system, we consider the reac- 
tion terms for the case where particles diffuse with diffusivity D and annihilate 
on collision. We shall see that a consequence of annihilation on collision is that 
g{Q) = for all t > 0, where g{y) is the correlation function defined in (|l7|). 
We derive an exact relation between ^'(0) and the rate of collisions between 
particles. 

Let s{y, At) be the probability that two particles, with initial separation y, 
collide before At. Then the reaction kernel K{y), defined in (^2|), is given by 



n§i 



If both particles diffuse with diffusivity D the: 

where we assume L ^ (DAt) i . To calculate the frequency of collisions between 
particles, we consider a time interval t,t + At. Given the density p(t) and the 
correlation function g(y, t) defined in ([l^), we can imagine following the paths of 
all the particles from time t to time t + At without removing those that collide. 
Then, the probability that a particle chosen at random undergoes a collision 
between time t and time t + At \s P{t, At), where 

P(t, At) = 2p{t) [ dy siy, At)g{y, t). (24) 
Jo 

The expression ( p^ overestimates the number of collisions in the system with 
annihilation on collision due to the possibility that the same particle undergoes 
two (or more) collisions in the interval i, t + At. However, this latter probability 
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is proportional to (At)^ as At 0, and so ( p^ is valid for our system in the 
limit At^O. 

Next, consider the dynamics of the system as a whole. The mean number 
of distinct collisions between time t and time t + At is given by Lp{t)P{t, At). 
We can, therefore, write 

/ dys{y,At)giy,t) 
Jo 



dy s{y, At) (^5(0, t) + 2/<?'(0+, t) + -y^g"iO+, t) + 

-^(8DAt)ig(0, t) + -(8DAt)g'(0+, t) + 0(DAt)i 
Vtt 4 

(25) 

The number of collisions between time t and t + At is proportional to At if 

5(0,i) = 0. (26) 
Because the number of nucleation events between time t and time t + At is pro- 



portional to 2LAt, the condition (26) is necessary if there is to be a steady-state 
balance between nucleation and annihilation. More generally, it is necessary if 
p(t) is to obey a differential equation. It is, of course, possible to construct ini- 
tial conditions that do not satisfy (p6|): a random distribution of particles, for 
example. Then the number of annihilation events will initially be m-QpQt:imial 



to 1 2 ; this period of rapid annihilation creates a "depletion zone '5@^^ in 
g{y,t), which thereafter satisfies (p6|). That g{y,t) < 1 for j/ — > implies an ef- 
fective repulsion: particles are less likely to be found close to a reference particle 
than a large distance from it. 

Using (^5|) and (^^ gives exact expressions for the evolution of the density: 

p{t) = -4:Dp'^{t)g'{0+,t) + 2T, paired; 

(27) 

p{t) = -ADp^{t)g'{{)+,t) + Q, unpaired. 

In particular, we have the following relationship between the steady state density 
and the derivative of the correlation function. Let po a-nd (7o(y) denote the steady 
state density and correlation function. Then 

T^2Dplg'^{Q+) paired; 

(28) 

Q = 4Dp§go(0+) unpaired. 
The reaction kernel K{y) is a singular function 

^ dy K{y)g{y,t)= lim^^^J^ dy s{y, At)g{y,t) ^ 2Dg'{0+ ,t). (29) 
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We have assumed that (M) holds. The derivative of g{y) is one-sided 



9'{0+,t)^ lim (30) 

a^0+ a 



because g{y) is only defined for y > 0. In other words 



Kiy) = 2D lim (31) 



It is interesting to, compare ( |3l| ) with the form of K{y) used, for example, by 
Lindenberg et alE3 

Kiy) = k5{y - a), (32) 

which introduced a reaction radius a and a rate coefficient k. There it was 
assumed that these constants are connected to the diffusivity via the Smolu- 
chowski relation k = 2D /a and that in the limit a — > 0, fc — > 00 their product 
remains finite. In the form used here, by contrast, we are able to explicitly take 
the limit of zero reaction radius: a — > 0. 
Similarly to (|§) 

y ^ dz if (|z|)F3(y, z; t) = ADFl,{y, 0+; i), (33) 



where 



and 



F!,{y,Q+-t) = lima-'F3{y,a;t), (34) 



L 

dz K{\z - y|)F3(y, z; t) = ADF^^iy, y+;t). (35) 

-L 

A similar expression was derived for a discrete, coagulation model without nu- 
cleation by Lin, Doering, and ben-Avraham.E3 Since F3{y,0;t) = Fz{y,y]t), 
Eqs. and ( |2l| ) simplify to the pair of equations 

p{t) ^ -ADp\t)g'{Q+,t) + 2T, (36) 

d d'^ 

Ql9{y,t) = 2D—g{y,t)-K(y)g{y,t)- 8Dp{t)F^{y,0+;t) 

+ 4rp-\t) + p-\t)q{y) - 2p-\t)g{y,t)p{t)- (37) 

Annihilation on collision is described by the terms involving the reaction kernel 
Kiy). 
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2.3 Truncation of the hierarchy 



We have obtained exact expressions for the evolution of the density. However, 
to obtain a closed set of equations, we truncate the hierarchy of distribution 
functions via an approximation. Various methods ]jiaKe|-heen used to break 
hierarchies resulting from reaction-diffusion systems .H'E3' Ell We shall restrict 
ourselves to the simplest. In the hierarchy that begins with ( ^6| ) and (|3^) we 
make the ansatz 



F^{y,0+;t)^g{y,t)g'{0+,t). (38) 

This choice, which would be exact if successive interparticle spacings were inde- 
pendentj^j is not per se the most compelling, but it has been shown to produce 
excellent results (when compared with siiruLlalions) for batch reactions and in 
the steady state with unpaired nucleation.EiEJ In Sec. || we shall compare the 
steady-state density obtained with this closure to the exact result. 



With the approximation (|38|) we find -4i:)p2(t)i^^(y,0+) = g{y,t){p{t)-2T), 
and so (|3^) and (|37| ) reduce to the following closed set of equations, linear in 

p{t) = ~ADp^{t)g'{Q+,t) + 2T, (39) 



^t9{y.t) = 2D^^g{y,t)~2Kiy)giy,t) + ^^il~giy,t)) + ^^, 



(40) 



plus the condition (26). In the case of paired nucleation, q{y) ~ rS{y — b). In 
the case of unpaired nucleation in the (thermodynamic) limit, L ^ p{t)^^ , ( ^o|) 
reduces to 



d 



d 



2Q 



^Q{y,t) = 2D—g{y,t)~2Kiy)giy,t) + ^il-g{y,t)) 



dt 



dy 



Pit) 



(41) 



with, as before, Q = 2r. Note that no further simplifications can be obtained 
by assuming low density. In particular, a low density expansion cannot be used 
to justify the truncation (38). 



2.4 Steady states 

The density and correlation function in the steady state, pq and go{y), are found 
by setting to zero the time derivatives on the left-hand side of (|3^) and of ( pO[ ) 
or (pl|). There thus results a second-order equation for go{y), with the two 
relations ( p^ ) and (p8|). 
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Figure 2: Correlation function for unpaired nucleation. Numerical results are 
compared with the formula (^3|) obtained by truncating the hierarchy (solid 
line). Q = 1.0 and D = 0.5. 



For unpaired nucleation one findsEil 

- ^ (42) 



V16L' 

9l-{y) = i-e-'(Q/'^^'"y, (43) 
where we have introduced the superscript "t" to indicate that the result is 



obtained from the truncation (38) and "u" to denote "unpaired nucleation." In 
Fig. ^ we compare numerical results for the correlation function with (^) . 
For paired nucleation the steady-state equation for the correlation function 

is 

= 2D^go(.y)^2Kiy)goiy) + i-[l-goiy)] + ^5{y-b). (44) 
dy^ Po Pa 

The solution of (|4^) is derived in Appendix A. When b is sufficiently large 
the results are equivalent to (^3|) with (|4^ ). Of interest here is the opposite 
situation: e — > 0, with e defined in Eq. (111). The separation b in the latter case 
is much smaller than the length scale defined by the inverse density, and the 
steady-state density is given by 
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Figure 3: Correlation function for paired nucleation. Numerical results are 
shown as dots and the approximation ( p6| ) as a solid line, b = 0.05, F = 64.0 
and D = 2.0. 



The correlation function in the same limit is 



y 



0<y<b- 
y>b. 



(46) 



Corrections to (^6|) are proportional to e * . In Fig. || this correlation function is 
compared with numerical results. 



2.5 Relaxation to the steady state 

In order to study the relaxation to the steady state, we decompose the functions 
p{t) and g{y,t) as follows: 



p{t) = po + Spit) 



(47) 



9o{y) + 5g{y.,t) , 



(48) 



with po a-iid go{y) the steady-state density and the steady-state correlation func- 
tion, respectively. This decomposition is valid for both unpaired and paired 
nucleation. Assuming that we are close to the steady state, we can obtain lin- 
earized equations for the deviations 6p and 5g from their steady-state values. 
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For paired nuclcation 
d 



Q^Sp{t) = -4Dpldg'{0+,t)-8Dpog'o{Q+)Sp{t) , (49) 



-Mv^i) = 2D^Sgiy,t)~2K{y)6g{y,t) 



4r 



Sg{y,t)+'—^6p{t) 

PO 



2r 

—S{y - b)Sp{t) 



Po 

(50) 



For unpaired nuclcation the last term in (50) is absent. 

Formal solution of these coupled linear equations is presented in Appendix B. 
Explicit solution for all times is in fact possible for the unpaired nuclcation 
case (and presented in the appendix). Ultimately we are interested in their 
asymptotic relaxation behavior. If the relaxation processes each involve a single 
exponential decay 

Sp{t) *^ Ae-"\ (51) 

Sg'{0+,t) Be-^*, (52) 

then the density and correlation function decay on the same time scale, i.e., 
(3 = a. For unpaired nuclcation we find from the exact result ( B.22| ) that the 
asymptotic decay is indeed exponential, with 

a" = (5 + %/5)(OT2)1/3^ 7.236...(OT2)1/3 

^ ^ '^(I?Q2)1/3 ^ 4.558 . . . {DQ'^Y'^. (53) 



41/3 

In the case of paired nuclcation, if we assume exponential decay we find for 
the inverse time scale 



1/2 



) • (54) 



Howeiier, on the basis of the exact results reported below (and also in parallel 
workEj) there is reason to suspect that the decay may not be purely exponential 
in the paired nuclcation case. 



3 Exact results 

In this section we derive exact expressions for the density of particles, using a 
function that obeys a linear partial differential equation. The function iSjSimilar 
in interpretation to the pair-pair correlation function in the Ising modelE3 The 
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methodology is also similar to. that used to obtain exact results for models of 
diffusion-limited coagulation£3 Here, we obtain explicit exact expressions for 
the density of particles, in steady state and nonsteady state, for paired and 
unpaired nucleation. Previous exact results for diffusion-lipiited reaction with 
annihilation have been limited to the case of no nucleationHj 
Let the function r{x, t) be defined as follows: 

r(x, t) = {probability that the number of particles (55) 
between and x at time t is even}. 

Note that, by translational invariance, we can replace the interval (0,x) by 
(X, X + x) for any X. The value of r{x, t) changes due to diffusion of particles 
in or out of of the region (0, x), and due to nucleation of a single particle in the 
region. 

In Appendix C we derive equations for the space and time derivatives of 
r(a;,t). At any time t the density p{t) is given by 

P(*) = -^Kx-^)L=o+- (56) 

To describe the time evolution of r{x,t), we distinguish between the cases of 
unpaired and paired nucleation using the superscripts u and p. 
In the case of unpaired nucleation, r^{x,t) satisfies 



cP" 

— r"(x,i) = 2i?^r"(a;,i)-xQr"(a;,i)+xQ(l-r"(a;,i)) 

= 2i?— r"(a;,i) +xQ(l-2r"(x,t)), (57) 



with the boundary conditions 



r"(0,i) = l and lim i) = ^, i > 0. (58) 



In the case of paired nucleation r^{x,t) satisfies 
d 



52 

2D—rP(x,t)+2xT(l-2rP(x,t)) x < b; 

ox^ 



2D—r^{x,t) + 2bT{1^2r^[x,t)) x > b, 

OX'' 



(59) 



with the boundary conditions 



rP(0,i) = l and lim rP(x, i) = i, t>0. (60) 
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Figure 4: The function rQ{x): numerical and exact results for unpaired nucle- 
ation. The sohd line is Eq. (|2|). Q = 1.0 and D = 0.5. 

3.1 Steady state: Unpaired nucleation 

The steady-state solution of ( ^7|) will be denoted by rQ{x). It satisfies 

2D^r^{x) + xQil - 2r^ix)) = 0. (61) 

The solution is0EI 

and is shown in Fig. ^. Thus the exact steady-state density for unpaired nucle- 
ation is 

Note that the exact result for the steady-state density is 0.9186 of the density 
( ^2|) predicted from the truncated hierarchy. For comparison, in a discrete model 
with nucleation rate R, where collision of particles prodiices coagulation rather 
than annihilation, the steady-state density is given byo'CJ 

„ I f |Ai'(0)| , , . , 
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Figure 5: The function Tq (x), measured at late times in a numerical simulation 
with paired nucleation (dots). The solid line is rQ{x) as given in Eq. (|65|). 
6 = 0.2, r = 0.25 and D = 0.5. 



3.2 Steady state: Paired nucleation 

The steady-state solution of (|9|) is 

i[ciAi((l^)^.)+c,Bi((f)^x) 



rP(x) 



x<b; 



X > b 



(65) 



and is shown in Fig. |5[ We have used the second of the boundary conditions 
( |60| ) to rule out increasing exponential solutions for x > b. The constants Ci, 
C2, and C3 are fixed by requiring rg(0) = and imposing continuity of rQ(x) and 
£rP(a;) at x ^ b. 

The density /9(t) is given by (|5^). In the steady state 



1 /2r\3 



(ciAi'(0) + C2Bi'(0)) 



4D 



^ |Ai'(0)| /Bi'(£) + V3Ai'(£)+£(Bi(£) + V3Ai(£)) 



Ai(0) VBi'(£) - V3Ai'(£) + £(Bi(£) - y3Ai(£)) 



(66) 



in terms of the dimensionless quantity e defined in Eq. (^. The function ( p6[ ) 
is plotted in Fig. ^. 
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2 3 

£ 



Figure 6: Exact steady-state density versus the dimensionless parameter e. 

In the limit e ^ 0, Bi(e) ^J'iAi{£) and Bi'(£) -y3Ai'(e), so 
p / r \^ / 1 1 



£2 £ 

AD) V 2 



2D 



(l + 0(£*)). 



For £ — > oo 



1 /2r\ ^ |Ai'(0)| 



2 V / Ai(0) ' 



(67) 



(68) 



and we regain the resuh (63) 



3.3 Time-dependent statistics: Unpaired nucleation 

Let us introduce 

/i^(a;,t) =r"(a;,t)-rS(a;). (69) 

Then, h^{x,t) satisfies 



o a2 

— h^{x, t) = 2D—h^{x, t) - 2Qx h^{x, t), 
at ax'^ 

with the boundary conditions 

/i"(0,t) = and lim /i"(a;,t) = 0, 



(70) 



(71) 
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for all t. We can expand the general solution as follows: 

oo 

h^{x,t) = J2c^K{^)e-^'', (72) 



i=l 



where the eigenfunctions hf{x) satisfy 

^Kix) - ^xKtix) = -\ih^{x), (73) 
with the boundary conditions 

/iJ'(O) = and lim ^^{x) = 0. (74) 

The eigenfunctions h^(x) are thus given by 

/^?W = iV.Ai((|)^x-A)j. (75) 
The eigenfunctions are normahzed by choosing 

^ = f ((§) ' - - if) ' II 

The eigenvalues Aj are related to the zeros of the Airy function (all on the 
negative real axis) 

Xi = -{8DQ^)^ai, (77) 

where is the ith zero counting away from 0. Relaxation towards the steady 
state is determined for late times by the smallest eigenvalue 

Xi = -{8DQ^)hi. (78) 

An explicit analytical solution for the density as a function of time is obtained 
once the constants c, are determined from the initial condition r^{x, 0) 



/>oo 

Ci= dx/i,"(x)(r"(x,0)-ro"(x)). (79) 
Jo 



Thus, 



K.(f)*,o)-^(i + A||) 



OO 

dz Ai{z + Oi) 

. g ^^^^^^^^^ 



^Ai{&^x-^)]e-^'\ (80) 
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i 


ai 


Ai'(a,) 


/-dzAi^(z) 


dz Ai(z - ai){l - 


Ai(^)^ 

" Ai(0) 


1 


-2.338 


0.701 


0.492 


0.972 




2 


-4.088 


-0.803 


0.645 


1.002 




3 


-5.520 


0.865 


0.749 


0.996 




4 


-6.786 


-0.911 


0.829 


1.001 




5 


-7.944 


0.947 


0.897 


1.012 





Table 1: Quantities related to the eigenvalues and eigenfunctions for unpaired 
nucleation. 



and 

p{t) 



dz Ai(z -I- tti) 



(81) 



/-dzAi^(z) 



X Ai'(ai)e' 



-Ait 



(82) 



3.3.1 Zero initial density 

If /9(0) = 0, then 0) = 1 for aU x > and 



Thus, 



Ai(z^ 



1 /n\^ °° /o°°dzAi(z-a,) l-T^ 



2\D ^ r°°dzAi^(z) 

In Fig. 1^, the exact time evolution is compared with numerical results, ob- 
tained with L = 3 X 10^. The lower dotted line is obtained by plotting only the 
first term of the sum in ( p5[ ) , using the values from Table |l|. Explicitly, the first 
eigenvalue that determines the long-time approach to the steady state is 

Ai = 4.676 . . . {DQ^f^ = 7.4227 . . . {DV^f/^. (86) 
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Figure 7: Time evolution starting with no particles present. Unpaired nucleation 
with Q = 16, D = 1. Solid circles are numerical results. The solid line, almost 
invisible under the numerical results, is the exact evolution calculated from 
Eq. (^. The upper dotted line is the exact steady state and the lower dotted 
line is the first term in the sum (|85|). 



3.3.2 Random initial density 

An interesting case is provided by starting the system with the exact steady- 
state density p(0) = pJJ, but with a random initial distribution of particles. 
There is an initial period of rapid annihilation that reduces the density, followed 
by a slower relaxation back to the steady-state value. 

For a random initial distribution of particles with density p, the number of 
particles in (0, x) is a Poisson random variable with mean px. The function 
r^{x,Q) can be calculated as follows: 

r"(a;,0) — 7'[0 particles between and x] 

+V [2 particles between and x] + . . . 

= l{l + e-'P^). (87) 

Figure^ shows data from a numerical simulation, performed with L — 2 x 10^, 
along with the results of the calculation of the coefficients in (|8^), using (|87|). 
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numerical 
exact steady state 
exact evolution 
first eigenfunction 



0.15 



Figure 8: Time evolution starting from a random distribution of particles with 
the exact steady state density. Unpaired nuclcation with Q = 16, D = 1. Solid 
circles are numerical results. The solid line is the exact evolution, the upper 
dotted line is the exact steady state, and the lower dotted line is the most slowly 
decaying term in the sum (|82| ) with initial conditions (87). 



3.4 Time-dependent statistics: Paired nucleation 

Let us introduce 

Then hP{x,t) satisfies 

d 



dt 



hP{x,t) ^ ChP{x,t), 



where the operator £ is defined by 



2D-^f{x) - ^xTf{x) x<b; 
ax 

2D-^fix) - 4&r/(x) x>b. 
ax 



The boundary conditions on K^ix^t) are 

h^{0,t) — and hP{x,t) bounded as a; ^ oo, 

for t > 0. 

Let us introduce 



(88) 



(89) 



(90) 



(91) 



Oil 



AbV 



(92) 
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For ai > I, the eigenvalue equation 

has a continuous spectrum of solutions 

' ci(a,) [V3 Ai (e(f - a,)) - Bi (e(f - a,))] x<b; 



(93) 



C2(Q;,)sin {a, - l)i {x - b) 



+ C3(a»)cos (a, - l)5(2rfc)5(a;~5) 



(94) 



> b. 



When e is sufficiently large, there are also discrete eigenvalues at values of < 1 
satisfying 

Ai(-ea,)Bi'(e(l - a,)) - Bi(-ea,)Ai'(e(l - a,)) 

+ {e{l - ai))^ (^Ai(-£ai)Bi(e(l - a^)) - Bi(-£a,;)Ai(e(l - a^)) 
= 0. 

The eigenfunctions in this case are 

ci(a,) [V3Ai (£(|-aO) - Bi (£(| - a,))] 



(95) 



X < b: 



(96) 



C4(a,)exp[-(^)2(1 -a^)■^x\ 



> b. 



Thus, for all finite e, there is a continuous spectrum of eigenvalues with 
Xi > 46r. For £ smaller than £c these are the only eigenvalues. The critical 
value Sc satisfies 



so 



Ai(-£,)Bi'(0) - Bi(-£,)Ai'(0) = 0, 



1.98. 



(97) 



(98) 



Discrete eigenvalues appear for larger values of e (Fig. The unpaired limit 
is regained as e ^ oo, when —eai — > a^, so that Xi 2.338(32Z)r^) a . 

In Figs. |l^ and |ll| we compare the exact expressions for the steady-state 
density and exponents characterizing relaxation toward the steady-state density 
with numerical results. In each case the curved dotted line is obtained using 
the lowest exponent available, but the coefficient is obtained from a best fit. In 



the case depicted in Fig. 10, there is a discrete eigenvalue Xi < ibT; in the case 
depicted in Fig. there is only the continuum of eigenvalues Xi > AbT, and 
the relaxation process may not be purely exponential. 
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A 




Figure 9: Eigenvalues for paired nucleation. All values A > 46r are permitted. 
Discrete values Aj < 46r are also found for sufficiently large e. The dotted line 
is A = -ai4bT/s. 
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Figure 10: Density of particles versus time for paired nucleation starting from 
zero density. The parameters are F = 16, D = 0.5 6 = 1 (e = 4). The solid 
circles are numerical simulation results and the upper dotted line is the exact 
steady state. The lower dotted curve is po{l — 1.6 exp(— 2.346rt)). Note that 
the latter exponent is the lowest for e = 4 and is a discrete value below the 
continuous spectrum. 
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Figure 11: Density of particles versus time for paired nucleation starting from 
zero density. The parameters are F = 16, I? = 0.5 & = 0.25 (e ~ 1). The 
sohd circles are numerical simulation results, the upper dotted line is the exact 
steady state. The lower dotted curve is po(l — 0.5 exp(— 46Ft)). The exponent 
chosen for the fit is the lowest in the continuum. No discrete eigenvalues are 
found for e = 1. 

3.5 Time-dependent statistics: No nucleation 

In the absence of nucleation, r"(a;,t) satisfies the heat equation 



with the boundary conditions 



r"(0,t) = 1 and lim r"(a;,t) 



(99) 



(100) 



for alH > 0. The solution of (|99| ) is given byS 

r-^{x,t) = l + iSirDt)-^ i^J dy (1 - r"(j/, 0))e-("-^''/«^* 

+ / dy(r"(-2/,0)-l)e"("-«)'/8^*j. (101) 

J — oo 

If the initial distribution of particles is random with density p(0), then 

r"(x,0) = i (l-e"2p(")-) . (102) 
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Now, using (pq), we derive the density of particles as a function of time for 
random initial conditions 

Pit) = -{8nDt)-^(^J^ aye~M0)tJL.^~yVsDt 







^ 8Dt 



1 f°° 

2p(0)(87ri:>i)~5 / d2/e-2K0)ye 

^0 



-2p{0)y^-y'/SDt 



= p(0) exp {8Dtp^{0)) erfc (p(0)(8Dt)^) . (103) 
We thus reproduce the result (^ of Torney and McConnell.0 

4 Discussion 

In the case of unpaired nucleation, there is only one length scale, proportional 
to {D/Q)3 ^ and only one time scale, proportional to {DQ^)^3 , The relaxation 
time to equilibrium and the mean lifetime of a particle are nmportional to 
one another. This is made clear in the mesoscopic approach,L3 although in 
the approaches detailed here we do not make this explicit distinction. It is 
noteworthy that the reaction-diffusion approach yields a steady state that is 
within 9% of the correct one and a relaxation rate that differs from the exact 
result by only 2.3%. 

Although the case of paired nucleation is more complicated in terms of the 
time scales associated with its dynamics, its steady-state distribution of particles 
is closer to a classical equilibrium random distribution than the corresponding 
distribution produced by unpaired nucleation. Note that the truncated hier- 
archy approach in this case leads to the exact steady state density. The un- 
derlying reason is that the dynamics produced by paired nucleation is close to 
time-reversal invariant. For comparison, an ensemble of noninteracting diffusing 
particles has a two-point function g{y) identically equal to iH We can imagine 
producing a spacetime diagram such as shown in Fig. |^ from a diagram associ- 
ated with noninteracting particles in two steps. First, when two particles collide, 
move them to a different, randomly chosen part of the system. Second, separate 
them by a distance b. The first step does not affect the correlation function or 
time-reversal invariance. The second step directly changes the correlation func- 
tion for separations smaller than b. Thus, for diffusion-limited annihilation with 
paired nucleation, as the parameter for e that measures the distance between 
newly nucleated pairs tends to 0, the two-point function in the steady state is 
appreciably different from 1 only in a region whose width is proportional to b. 

Results obtained for the steady-state density and relaxation rate are summa- 
rized in Table 0. For the evaluation of time scales, the case of paired nucleation 
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Unpaired Nucleation 


Paired Nucleation 




b OO 


5^0 




Po 


Relaxation rate 


Po 


Relaxation rate 


Mesoscopic 


cx iQ/Df' 


oc(Z?Q2)i/3 


cx {bV/Df'^ 


cx [b/TDf^, bV 


Hierarchy 


{Q/l&Df'^ 


0.219/ (DQ2)i/3 


{bT/2Df''^ 


{b/mDYf^ 


Exact 


0.9186 {Q/lQDf/^ 


0.2138/ [DQ^y ' 


{bTl2Df'^ 


46r 



Table 2: Summary of results for steady state densities and relaxation rates 
obtained by various methods. 



is more complicated. The different approaches indicate the occurrence of multi- 
ple length and time scales. The mesoscopic approachcHl leads to a characteristic 
time for approach to equilibrium proportional to (bV)^^ , and a distinct mean 
lifetime of a particle proportional to {b/ DT)"^ . These two time scales were iden- 
tified as corresponding, respectively, to recombination (two particles created a 
distance b apart collide and annihilate) and to nonrecombinant annihilation (col- 
lision between two particles nucleated at different times). Unpaired annihilation 
is less frequent than paired|-apnihilation, but both time scales are important in 
the dynamics of the system.E2l The exact approach leads to the former time scale 
as an upper bound of a continuum of scales. The truncated reaction-diffusion 
hierarchy leads to the latter time scale under the assumption of exponential 
decay, which may not be valid. Understanding the time scales in the case of 
paired nucleation requires further research. 

The theoretical approach based on truncation of a hierarchy of distribution 
functions thus permits the calculation of steady-state densities and correlation 
functions that are in fair agreement with simulations. It gives the exact result 
for the steady-state density in the limit £ — > 0, where the statistical distribution 
of particles is close to random. As pointed out by van Kampen,0 an approach 
using a truncation can be made systematic if it is based on an expansion in 
a small parameter. However, his suggestion that the small parameter be the 
density of particles is not applicable to the case of unpaired nucleation of point 
particles because there is no other quantity with the dimension of length, i.e., 
nothing for the density to be small compared to. The exact approach based on 
the function r{x,t) sidesteps these difficulties by providing a direct method to 
calculate the density of particles. For any value of £, the density in the steady 
state and its time-dependent statistics can be exactly calculated. However, the 
method has not yet been extended to exact calculation of the full distribution 
of interparticle distances. 
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APPENDIX A: Steady state solution of truncated hierarchy 

In this Appendix wc find the steady state solutions of ( |40| ) and (^l]). 

Unpaired nucleation We solve with the left-hand side set to zero. 
Fourier transforming the quasilinear equation according to 

dy go{y)cos——, (A.l) 



^ ' Jo 



leads for n 7^ to 



2r 



^ ^" (A.2) 



SDn^n^ , 4r 



where Eq. ( p8| ) has been used. Normalization sets go = i^"'^^^. Fourier inversion 
according to 

, , 1 ^ 2Trny 

ayyi^JTH 1^ 5nC0S— — , (A.3) 



n— — oo 



can be done by separating out the n = Q contribution explicitly and changing 
the resulting sum to an integral (valid as L oo) 

go[y) = 1 2 / KjT- 



T^Po Jo q2 



Dpo 



1 / 2r 



1/2 



4 \Dpl} ^ > 



The requirement ( |26|) must also hold in the steady state. Together with (A.4) 
this leads to Eqs. @ and (|3|). 

Paired nucleation We now solve Eq. (p4[), where the steady-state density po 
is related to the derivative of go [y] as j/ ^ 

r = 2Dg'^{0+)pl. (A.5) 

The (5-function contribution in the last term of leads to a discontinuity in 
the derivative g'^^y) a,t y — b. This leads to the search for a solution of the form 

l + (5-l)e-(2r/^''«)'''?' 0<y<b; 
go{y) = { (A.6) 
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The constant 5' is determined from the condition (A. 5) 

1 

( r 



1 - 25 = 



(A.7) 



Now, the constant P is determined by enforcing continuity of the solution .go(y) 



P = S{1~ e2(2r/i?Po)^/^b) _ i_ 
The discontinuity in the derivative gQ^y) at y = b is 

r 



j'o{b+)-9'o{b-) = - 



2Dpf, 



2 ■ 



(A.; 



(A.9) 



Using ( [A.6| ) to evaluate the left-hand side of (A.9) and rearranging gives an 
implicit expression for the steady-state density 



1/2 



or, rearranging again 
where u is defined by 



8L> 



1/2 



_L = _L _ 3/2\ (2r/Dpo)^/^. 



/863r\" 



(A.IO) 



(A.ll) 



(A.12) 



While we cannot invert (A.ll) explicitly for the steady-state density, we 
can examine the limits in the dimensionless parameter a. In the limit cr ^ 1, 
corresponding to large initial separation, we find 



3/2 



(A.13) 



so Po r/SL*, as obtained in ( ^ ) for the unpaired nucleation case. Morevover, 
S" ^ and g{y) 1 - e"^''"^ as in (||). 

The limit cr ^ corresponds to small initial separation. Then (A.ll) reduces 



to 



(A.14) 

Note that cr ^ corresponds to bp^ — > 0. Expanding in powers of ct, we find 

5" = 



Po 



bV 
2D 



1 



(l--a + 0(a2)), 



' ^ ' 0(ct) and P=icr + 0( 



16 



2 

(A.15) 
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The correlation function go can be expanded as 

f(l + ia)+0(a2) 0<y<6; 
9o{y) - { (A.16) 
l + \<Te-''y/^ + 0[a^) y>b. 

These results as u — > are reported in Eqs. (Esh and (Eq). 



33 



APPENDIX B: Relaxation to the Steady State 

Here we detail the calculation of approach to the steady state in the trun- 
cated hierarchy approach. 



Unpaired nucleation It is convenient to introduce the symbols 



(B.l) 



The linearized perturbation equations ( [49| ) and (50) in the unpaired case are 
d 



dt 



d_ 

dt 



5p{t) = -x[V(0+,t) + 8<5p(i)] , 



B.2) 



-Mv^t) = 2D-^6giy,t)~2Kiy)Sgiy,t)-8xSgiy,t) 
-I67 [l~g^{y)]6p{t), 



(B.3) 



where the specific form (|4^ ) has been implemented. As initial conditions we 
choose 



6g{y,t^ 0)^0, 



(B.4) 



and an arbitrary Sp{0). Note that having implemented the condition ( p^ ) on 
the steady-state solution implies that 6gn=o{t) = 0. 
The solution of Eq. ^3) is 



(B.5) 



Sp{t) = e-'^^'SpiO) - X / dre-«'^(*-"Ug'(0+,T) 



Transforming (|B.3| ) according to (Al) gives, for n 7^ 

^^^9n{t) - - ^^^2" ^9n{t) - SxSgnit) - ^6g'{0+,t) + m-fgnSpit). (B.6) 
We can formally solve this equation as well, to obtain 



dre 







^5g'{Q+,T) + lQr9n5p{T) 



(B.7) 



and hence its inverse Fourier transform 







dr /C(y, t - r) V(0+, r) + 8 / dr ^(y, t - t)5p{t). (B. 
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We have introduced the following functions (and taken the limit L — > oo): 

I +00 



K.{y,t) 



8D 



E 



2TTny s-dtt 
cos e L 

n— — 00 L 



8D 



(27rL'i)i/2 
27 



3 2Dt ^ 



8xt 



/ +00 

2Tmy _ aoj, 
cos — - — e i 

n— — 00 Jj 



E 



5n 



(B.9) 



872 _ 



8xt 



dq e 



-2Dq^t_ 



8x + 2Dq^ 



cosqy, (B.IO) 



where the prime on the sums indicates omission of the n = term and we have 
used Eq. ( |Al2| ). 

Since the unknowns appear on both sides of ( |B.5| ) and ( |B.8| ), these are only 
formal solutions. To proceed, we Laplace transform them (indicated by a tilde, 
Cf{t) = /(s), and solve the resulting set self-consistently. The limit ( pO[) must 
be handled carefully and not implemented prematurely. We find 



(B.ll) 



IC{y,.s)S~g'iO+,s) + 8giy,.s)Spis) 



(B.12) 



From these two equations we can obtain an expression for 5g'{0^ , s) as follows. 
First, set y = in ( p3.12 ). Since g{0,t) — and since 5o(0) = 0, it follows from 
Eq. @) that Sg{0, i) = for aU t and therefore 6g{0, s) = for all s. Thus, we 
find 



= -/C(0, s)dg'iO+, s) + 8g{0, s)Sp{s). 



Subsitution of (B.ll) into this result immediately leads to 

(w + 1)^(0, s)+^(0,s)' 



(B.13) 



(B.14) 



where we have set u = s/8x- This, together with Eqs. ( |B.ll|) and (|B.12| ), 
constitutes a complete solution of the Laplace transform of the problem. The 
indicated transforms can be calculated explicitly 



^(0,s) 



2{D/V) 



1/3 



Vim 

4u 



1 

\/-it + 1 



(B.15) 
(B.16) 
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and one readily obtains 

Spis) = f - , . (B.17) 

X \i - VTTTi + 8u{i + u) J 

It is possible to Laplace invert these expressions exactly .El In particular, Eq. ( B.17 ) 
can be rewritten as 

m.) = '-^( +^ +^ V (B.18) 

X \Vu + 1 - yi Vu + 1 - 2/2 Vu + l-y^J 



where 



and 



yi = 

The inversion 



Ai = A2 = A3 = — (B.19) 

40 40 4 ^ ' 



i(V5-l) y2 = -i(V5 + l) 2/3 = -^. (B.20) 



£-1 ( , " ) = e-^* ( + Be^'*erfc(-BVi) I (B.21) 

can then be applied to obtain 

Sp{t) - 8(5p(0)e-«^* ^ Ay,e''^y"'elic{-y^^/8^t). (B.22) 

It is noteworthy that this solution is exact within the truncation approximations 
of the model, that is, it represents the full time-dependent solution for the model. 

Asymptotic analysis of the exact result yields pure exponential decay as 
indicated in Eqs. (51) and (p^), with 

Q = (5 + %/5)x 7.236.. .X- (B.23) 

The proportionality of 5g and Sp clearly leads to the same decay rate for Sg{y, t) 
as for 5p{t). 

Paired nucleation Here it is convenient to introduce the symbols 

4r 

n = —, e = vb. (B.24) 



The linearized perturbation equations ( }49| ) and ( |50[ ) in the paired case are 

^^5p[t) = -4DplSg'{0+,t)-nSp{t), (B.25) 

-<55(y,i) = 2D—6giy,t)~2K{y)dgiy,t)-n6g{y,t) 

4r 2r 

[1 - g^oiy)] Sp{t) - —S{y - b)6p{t). (B.26) 

Po Po 
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With the same initial condition as in the unpaired case, the solution of 

(B.27) 



Eq. (B.25) is formally given by 



Transforming ( B.26| ) according to (A.l) gives for n 



put) 



SD 



9n6p{t) - ^ , COS 



5p{t). 



(B.28) 



We can formally solve this equation as well, to obtain 



dre 



x{-|^W,-) + 



n 



4T (2Tinb 
■ cos 



5p{T) 



and hence its inverse Fourier transform (notice that (5gn=o(0 = fo'" ^) 



dTlC{y,t-T)5g'{Q+,T)+ / dTg{y,t~T)5p{T) 



4r r 
p3fTj2 dTH{y,t-T)6p{T). 



(B.29) 



We have introduced the following functions 

/ +00 



/C(y,t) 



8D 
IT' 



-fit 



cos ■ 

= ~oo L 



2'Kny soTT^ n^ ^ 



8D 



e e 2Dt , 



Giy,t) 



n{y,t) 



PoL'/^ 
1 



27rny sp^l 
cos — - — e ^ 

n— — cxD 



E 



(B.30) 



(B.31) 



ot 



, +00 

2Trny 2'Knb _mj!^^t 

cos — - — cos — - — e . (B.32) 

n= — oo Li L 



E 



We Laplace transform ( p3.27[ ) and (B.2S) and solve the resulting set self- 
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consistently, to obtain 



59{y^s) = -IC{y,s)6g'{0+,s) +g{y,s)6p{s) 
from which in turn we find 



4r 



-IC{y,s)Sg'{0+,s) 



4r 



S{y^s)--^^n{y,s) 



6p{0)-ADpl- 



(B.33) 



We are interested in obtaining 5^'(0+, s). This can readily be done by eval- 
uating the previous equation at y = 0, to obtain (notice that we have chosen 
6g{y = 0,t) = for all t) 



= -iC{0,s)Sg'{0+,s) + 



1 



-Sp{0) - 4Dpl 



S~g'{0+,s) 



_s + n '^^ ' s + Ji 

The previously introduced functions evaluated at the origin become 



/C(0,t) = 

a(o,t) = 

W(0,t) = 



-Of 



+ 00 



Qt 



1 



n=—oo 



+ 00 



2'Knb _SDlllnlf 

cos — - — e 



(B.34) 



(B.35) 



(B.36) 



(B.37) 



We are interested in obtaining the Laplace transform of these functions. It is 
easy to obtain 



2(D/r)i/3 



/C(0, s) 

n{0,t) = i^e-"* / dge-2^9 * cosg6, 

Jo 



^s/n + 1 ' 

LV2 



(B.38) 

(B.39) 
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so that 



4r 



2r 



__^-nt^-b'/(sDt) 



4r 



172^(0, 



2r 1 -b,/i±ir 

—— p "V 2D 



(B.40) 



(B.41) 



The Fourier components for the two-point correlation function associated 
with Eq. ( A. 6 ) are 



2 v/2r/(ZJpo) 

rl/2 2r I 47r^7t^ 
Z5po + 



(25 - 1) - 25eV2r/(i?Po) 



2'Knh 



cos ■ 



(B.42) 



This form in turn leads to the Laplace transform of Q(Q,t): 



1 




COS qb 



2r/{Dpo)){s + n + 2Dq^)' 



(B.43) 



It is now a straightforward matter to 1) collect the various Laplace transform 
expr ession s to solve for Sg'{Q^, s) using Eq. ( B.34 ), 2) substitute this result into 
Eq. ( B.ll ), and explore the poles of the denominator of the resulting 6j5{s) in 
the limit e — *■ 0. The procedure is tedious but leads to the inverse time scale 
(p4|). The proportionality of Sg and 6p clearly leads to the same decay rate for 
6g{y,t) as for Sp{t). We do note that it is not clear from this procedure that 
the relaxation process is actually exponential in time. If it is exponential (and 
there is reason to question this from the results of the exact and mesoscopic 
procedures) , then it is necessary to perform the inverse Laplace transform more 
carefully. This is possible, but beyond the scope of this paper. 
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X X + Ax 



Figure 12: The intervals (0,a;) and (x,x + Ax). 



APPENDIX C: The function r{x,t) and its derivatives 

We first calculate the derivative of r(a;, t) with respect to x by considering 
the intervals shown in Fig. 

• Let Vcoix, lS.x^t) be the probability that there is an even number of par- 
ticles in (0, x) and no particle in (x, x + Ax) at time t. 

• Let Vci{x, Ax,t) be the probability that there is an even number of par- 
ticles in (0, x) and one particle in (x, x -I- Ax) at time t. 

• Let Voi{x, Ax, t) be the probability that there is an odd number of parti- 
cles in (0, x) and one particle in (x, x + Ax) at time t. 

The function r{x,t) defined in Eq. ( |55| ) can be expressed in terms of these 
quantities as follows: 



r(x-fAx,t) = Vco{x,Ax,t)+Voi{x,Ax,t) + 0{Ax^) 

r{x,t) = Veo{x,Ax,t)+Vci{x,Ax,t) + 0{Ax^). (C.l) 

Thus 

r(x + Ax,t) ~ r[x,t) = Voi[x,Ax,t) - Vci[x, Ax,t) + 0[Ax^). (C.2) 
In particular, by choosing x = 0, 



r{Ax,t) - r{0,t) 



0-Voi{0,Ax,t) + O{Ax^) 



Axp{t) + 0{Ax^). 



(C.3) 



Thus, the density p{t) is given by 



Pit) 




a;=0+' 



(C.4) 



which proves Eq. (|5^. 
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• k- 

I \ \ 1 

X — Ax X X + Ax 



-H \ 1 

X — Ax X X + Ax 



Figure 13: Particle movements contributing to the change in r{x,t). The effect 
in each case depends on whether the number of particles in (0, x — Ax) is odd 
or even. 

Next, we consider the intervals {0,x — Ax), {x — Ax,x), and (x,x + Ax). 
Let Veoi{x, Ax, t) be the probability that there is an even number of particles in 
(0, X — Ax), no particle in {x — Ax, x), and one particle in {x, x + Ax) at time t. 
Let Vem{x,Ax,t), Vcwix, Ax,t), Vcii{x, Ax,t), Vom{x,Ax,t), Vooi{x, Ax,t), 
'^010(2;, Ax,t), and Vouix, Ax,t)he defined in the obvious way. The appropri- 
ate intervals are shown in Fig. |l3| . 

Because g{0,t) = for all t > 0, the probability that there are two particles 
in (x — Ax, X + Ax) is proportional to Ax^ as Ax — > 0. Thus, 

r{x + Ax, t) + r[x - Ax, t) - 2r{x, t) 

= 7^010(2;, Ax, t) - Vowix, Ax, t) 

-Vcoiix, Ax, t) + Vooi {x. Ax, t) + 0(Ax3). (C.5) 

We derive the contribution due to diffusion of particles to the partial dif- 
ferential equations (|5^) and (|5|) for the evolution of r(x, t) by considering the 
probability that a particle at x — Ax at time t diffuses out of the region (0, x) 
before time t -\- At, and the probability that a particle at x -t- Ax at time t 
diffuses out of the region (0, x) before time t -t- At. 

The probability that a particle, at x-l- Ax at time t, is in (0, x) at time t + At 
is given byl23'L3 



1 r°° 1 / Ax 

Q{Ax,At)^{ATTD) 5 / dx expi-x^ /ADAt) = -erfc 

J Ax 2 V(4£'At)2 

(C.6) 
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Let 



d 

Rcoi{x,Ax,t) = ——'Pcoi{x,Ax,t) 
oAx 



(C.7) 



As Ax —^ 0, this quantity is the probabiHty density for finding a particle at 
X + Ax at time t, given that the number of particles in (0, x) is even. Similarly, 
let 



d 



d 



Rew{x,Ax,t) = -^-j;^Vcio{x,Ax,t), Rooiix,Ax,t) = -^^Vooi{x, Ax,t), 



dAx 

and so on. 

The time derivative of r(a;, t) 



{CI 



d 1 
— t) = ^im^ At ^ + ^ 



(C.9) 



is found using Eq. (C.5), integrating over Ax, and taking the limit — > 

/•OO 

r{x,t + At)~r{x,t) = 2 dAx Q{Ax, At) {Rcio{x, Ax,t) ~ Roio{x, Ax,t)) 

Jo 

+2 / dAx Q{Ax,At){-R^oi{x,Ax,t) + Rooi{x,Ax,t)) 
Jo 

= 2 J dAxQiAx,At)-^[raio{x,Ax,t)-VMx,Ax,t) 



-Vcoi [x, Ax, t) + Vooi {x, Ax, t) 



2 / dAx Q{Ax, At)—— —^r{x, t)Ax^ + 0{A 
Jq a Ax \ax'^ 

L 



°° f Ax 

dAx erfc 



\{ADAt)'2 J \dx^ 



;r{x,t)Ax + 0{Ax'^) 



2D—r{x,t)At + 0{At^) 



(C.IO) 
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